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ABSTRACT 

39 galaxies are now known, from follow-up of faint IRAS sources and from 
submm observations of high redshift AGN, with far infrared luminosities > 
lO 13 ^©. 13 of these, which have been found in 60 or 850 /im surveys, form 
an important unbiased sub-sample. 12 have been found by comparison of 
60 /im surveys with quasar or radio-galaxy catalogues, or from infrared 
surveys with colour selection biased towards AGN, while a further 14 have 
been found through submm observations of known high redshift AGN. In 
this paper I argue, on the basis of detailed modelling of the spectral energy 
distributions of hyperluminous galaxies with accurate radiative transfer 
models, and from evidence of high gas-mass in several cases, that the bulk 
of the emission from these galaxies at rest-frame wavelengths > 50/im 
is due to star formation. Even after correction for the effects of lensing, 
hyperluminous galaxies with emission peaking at rest-frame wavelengths 
> 50/im are therefore undergoing star-formation at rates > 10 3 M Q yr _1 
and are strong candidates for being primeval galaxies, in the process of a 
major episode of star formation. 

Key words: infrared: galaxies - galaxies: evolution - star:formation - 
galaxies: starburst - cosmology: observations 



1 INTRODUCTION: HAVE WE 

DETECTED PRIMEVAL GALAXIES ? 

This paper is directed towards the question; Have we 
already detected primeval galaxies ? The characteris- 
tics of a primeval galaxy might be 

• high redshift (z > 1) 

• undergoing a major episode of star formation 
(to form 1O 11 M0 in < 10 9 yrs, we need a star for- 
mation rate M* > W 2 M Q yr~ 1 ) 

• high gas fraction, say M gas ~ 10 Mq 

• evidence of interactions, mergers, dynamical 
youth. 

First efforts to find such galaxies centred on 
spectroscopic searches for Lya-emitting galaxies (see 
eg the review by Djorgovski and Thompson 1992). 
While examples of such galaxies are now being found, 
these early surveys suggested that either star forma- 
tion must be a more protracted process, occurring in 
smaller bursts (as expected in many bottom-up sce- 
narios), or that dust extinction must play a large role. 



Steidel et al (1996) have shown that star-forming 
galaxies at z > 3 can be found through deep ground- 
based photometry in the U, G and R bands. The 
high redshift galaxies manifest themselves as U-band 
'dropouts' as the Lyman limit absorption is redshifted 
into the U band. Over 500 spectroscopically confirmed 
high redshift galaxies have now been found by this 
technique. Many have weak or non-existent Lyman a 
emission, which accounts for the lack of success of the 
spectroscopic surveys. The role of dust in these galax- 
ies has been discussed by Pettini et al (1997, 1998a, b), 
Meurer et al (1997, 1998), Calzetti (1998). Pettini et 
al (1997, 1998a,b), Dickinson (1998) and Steidel et al 
(1998). Even at redshift 3 it appears that dust ex- 
tinction may be appreciable. However star formation 
rates in these galaxies are not exceptional, typically 
1-30 Mo/yr. 

The first evidence for galaxies with very high 
rates of star formation came from infrared surveys. 
Starburst galaxies had been first identified by Weed- 
man (1975) from their ultraviolet excesses and char- 
acteristic emission-line spectra. Prior to the launch 
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of IRAS, balloon and airborne measurements had 
demonstrated that Lji r > L opt for several starburst 
galaxies (see review by Sanders and Mirabel 1996). 
Joseph et al (1984) proposed that interactions and 
mergers might play a role in triggering starbursts. 
One of the major discoveries of the IRAS mission 
was the existence of ultraluminous infrared galaxies, 
galaxies with Lf ir > lQ 12 h^ Lq (/i 50 = H o /50). The 
peculiar Seyfert 2 galaxy Arp 220 was recognised as 
having an exceptional far infrared luminosity early in 
the mission (Soifer et al 1984). 

The conversion from far infrared luminosity to 
star formation rate has been discussed by many au- 
thors (eg Scoville and Young 1983, Thronson and 
Telesco 1986, Rowan- Robinson et al 1997). Rowan- 
Robinson (1999) has given an updated estimate of 
how the star-formation rate can be derived from the 
far infrared luminosity, finding 

M*,au/lL 6 o/L e ] = 2.2 <)>/e xlCT 10 

where cj> takes account of the uncertainty in the 
IMF (= 1, for a standard Salpeter function) and e 
is the fraction of uv light absorbed by dust, esti- 
mated to by 2/3 for starburst galaxies (Calzetti 1998). 
We see that the star-formation rates in ultraluminous 
galaxies are > 10 2 M Q ?/r _1 . However the time-scale 
of luminous starbursts may be in the range 10 7 - 10 s 
yrs (Goldader et al 1997), so the total mass of stars 
formed in the episode may typically be only 10% of 
the mass of a galaxy. 

In this paper I discuss an even more extreme 
class of infrared galaxy, hyperluminous infrared galax- 
ies, which I define to be those with rest-frame in- 
frared (1-1000 fj,m) luminosities, Ltoi,ir, in excess 
of lO 13 ' 22 /i M) 2 L (=10 13 '°/i 65 2 L o )- For a galaxy with 
an M82-like starburst spectrum this corresponds to 
1/60 > 10 13 ft"5o 2 Lq, since the bolometric correc- 
tion at 60 fj,m is 1.63. Sanders and Mirabel (1996) 
have a slightly more stringent criterion, Ltoi,ir > 
lO 13 /if 5 2 I/0, but in practice they use an estimate of 
L bo i based on IRAS fluxes, which results in a demar- 
cation almost identical to that adopted here. While 
the emission at rest-frame wavelengths 3-30 /j,m in 
these galaxies is often due to an AGN dust torus 
(see below), I argue that their emission at rest- frame 
wavelengths > 50/im is primarily due to extreme star- 
bursts, implying star formation rates in excess of 1000 
M jyr. These then are excellent candidates for being 
primeval galaxies, galaxies undergoing a very major 
episode of star formation. 

A preliminary version of these arguments was 
given by Rowan- Robinson (1996). Granato et al 
(1996) modelled the seds of 4 hyperluminous galaxies 
(F10214, H1413, P09104 and F15307) in terms of an 
AGN dust torus model. Hughes et al (1997) argued 
that hyperluminous galaxies can not be thought of as 
primeval galaxies on the basis of their estimates of the 
gas mass and star formation rates in these galaxies. 
However I show below that their arguments are not 
compelling. Fabian et al (1994, 1998) argued from X- 



ray evidence that 09104+4109 and 15307+325 are ob- 
scured AGN and this interpretation is supported for 
these objects by the non-detection of CO and submm 
continuum radiation (Yun and Scoville 1999, Evans 
et al 1999). However the very severe upper limits on 
X-ray emission set by Wilman et al (1999) for sev- 
eral hyperluminous infrared galaies led the latter to 
conclude that the objects might be powered by star- 
bursts. McMahon et al (1999) interpret the submil- 
limctrc emission from high redshift quasars and other 
hyperluminous infrared galaxies as powered, at least 
in part, by the AGN rather than by star formation. 
On the other hand Frayer et al (1999a,b) favour a 
starburst interpretation for 14011+0252 and 02399- 
0136. I discuss these arguments further below and try 
to resolve the question of what fraction of the far in- 
frared luminosity is powered by a starburst or AGN. 
The use of model spectral energy distributions derived 
from accurate radiative transfer codes is a significant 
advance on some previous work. 

I assume throughout that H — 50, 0. o = 1. With 
lower Q , more galaxies, especially at higher redshifts, 
would satisfy the definition adopted. 



2 PROPERTIES OF ULTRALUMINOUS 
INFRARED GALAXIES 

Sanders et al (1988) discussed the properties of 10 
IRAS ultraluminous galaxies with 60 /im fluxes > 5 
Jy and concluded that (a) all were interacting, merg- 
ing or had peculiar morphologies, (b) all had AGN 
line spectra. On the other hand Leech et al (1989) 
found that only 2 of their sample of 6 ultraluminous 
IRAS galaxies had an AGN line spectrum. Leech et al 
(1994) found that 67 % of a much larger sample (42) 
of ultraluminous galaxies were interacting, merging or 
peculiar. Lawrence et al 1989 had found a much lower 
fraction amongst galaxies of high but less extreme in- 
frared luminosity. The incidence of interacting, merg- 
ing or peculiar galaxies by ir luminosity is summarised 
in Fig 1 of Rowan- Robinson (1991): the proportion 
of galaxies which are peculiar, interacting or merging 
increases steadily from 10-20% at low ir luminosities 
to > 80% for ultraluminous ir galaxies. The situa- 
tion on point (b) remains controversial, though, since 
Lawrence et al (1999) find only a fraction 21% of 81 
ultraluminous galaxies in the QDOT sample to have 
AGN spectra (but on the basis only of low-resolution 
spectra). Veilleux et al (1995) find for a smaller sam- 
ple (21 galaxies) that 33% of ultraluminous galax- 
ies are Seyfert 1 or 2, with an additional 29% hav- 
ing liner spectra, which they also classify as AGN. 
Veilleux et al (1999a) find that 24% of 77 galaxies 
with 10 12 < L ir < 10 12 ' 3 are Seyfert 1 or 2, increas- 
ing to 49% of 31 galaxies with 10 12 ' 3 < L ir < 10 12 ' 8 
(for H = 75). Veilleux et al (1999b) find, from near- 
ir imaging studies, no evidence that liners should 
be considered to be AGN. Sanders (1999) reports 
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that most nearby ultraluminous ir galaxies contain 
an AGN at some level. 

Rowan-Robinson and Crawford (1989) found 
that their standard starburst galaxy model gave an 
excellent fit to the far infrared spectrum of Mk 231, 
an archetypal ultraluminous ir galaxy. However their 
models for Arp 220 appeared to require a much 
higher optical depth in dust than the typical star- 
burst galaxy. Condon et al (1991) showed that the ra- 
dio properties of most ultraluminous ir galaxies were 
consistent with a starburst model and argued that 
many of these galaxies required an exceptionally high 
optical depth. This suggestion was confirmed by the 
detailed models of Rowan-Robinson and Efstathiou 
(1993) for the far infrared spectra of the Condon et 
al sample. 

Quasars and Seyfert galaxies, on the other hand, 
tend to show a characteristic mid infrared contin- 
uum, broadly flat in vS v from 3 — 30pim. This compo- 
nent was modelled by Rowan- Robinson and Crawford 
(1989) as dust in the narrow-line region of the AGN 
with a density distribution n(r) a r _1 . More realistic 
models of this component based on a toroidal geome- 
try are given by Pier and Krolik (1992), Granato and 
Danese (1994), Rowan-Robinson (1995), Efstathiou 
and Rowan- Robinson (1995). Rowan- Robinson (1995) 
suggested that most quasars contain both (far ir) star- 
bursts and (mid ir) components due to (toroidal) dust 
in the narrow line region. 

Rowan- Robinson and Crawford (1989) were able 
to fit the IRAS colours and spectral energy distribu- 
tions of galaxies detected in all 4 IRAS bands with 
a mixture of 3 components, emission from interstel- 
lar dust ('cirrus'), a starburst and an AGN dust torus. 
Recently Xu et al (1998) have shown that the same 3- 
component approach can be used to fit the ISO-SWS 
spectra of a large sample of galaxies. To accomodate 
the Condon et al (1991) and Rowan-Robinson and 
Efstathiou (1993) evidence for higher optical depth 
starbursts, Ferrigno et al (1999) have extended the 
Rowan-Robinson and Crawford (1989) analysis to in- 
clude a fourth component, an Arp220-like, high opti- 
cal depth starburst, for galaxies with log Leo > 12. 
Efstathiou et al (1999) have given improved radiative 
transfer models for starbursts as a function of the age 
of the starburst, for a range of initial dust optical 
depths. 

Sanders et al (1988) proposed, on the basis of 
spectroscopic arguments for a sample of 10 objects, 
that all ultraluminous infrared galaxies contain an 
AGN and that the far infrared emission is powered by 
this. Sanders et al (1989) proposed a specific model, 
in the context of a discussion of the infrared emission 
from PG quasars, that the far infrared emission comes 
from the outer parts of a warped disk surrounding 
the AGN. This is a difficult hypothesis to disprove, 
because if an arbitrary density distribution of dust is 
allowed at large distances from the AGN, then any 
far infrared spectral energy distribution could in fact 



be generated. In this paper I consider whether the 
AGN dust torus model of Rowan- Robinson (1995) can 
be extended naturally to explain the far infrared and 
submilllimetre emission from hyperluminous infrared 
galaxies, but conclude that in many cases this does 
not give a satisfactory explanation. I also place con- 
siderable weight on whether molecular gas is detected 
in the objects via CO lines. 

Rigopoulou et al (1996) observed a sample of 
ultraluminous infrared galaxies from the IRAS 5 
Jy sample at submillimetre wavelengths, with the 
JCMT, and at X-ray wavelengths, with ROSAT. They 
found that most of the far infrared and submillime- 
tre spectra were fitted well with the starburst model 
of Rowan- Robinson and Efstathiou (1993). The ratio 
of bolometric luminosities at 1 keV and 60 /im lie in 
the range 10~ 5 — 10 -4 and are consistent with a star- 
burst interpretation of the X-ray emission in almost 
all cases. Even more conclusively, Lutz et al (1996) 
and Genzel et al (1998) have used ISO-SWS spec- 
troscopy to show that the majority of ultraluminous 
ir galaxies are powered by a starburst rather than an 
AGN. 



3 HYPERLUMINOUS INFRARED 
GALAXIES 

In 1988 Kleinmann et al identified P09104+4109 with 
a z = 0.44 galaxy, implying a total far infrared lumi- 
nosity of 1.5xl0 13 , a factor 3 higher than any other 
ultraluminous galaxy seen to that date. In 1991, as 
part of a programme of systematic identification and 
spectroscopy of a sample of 3400 IRAS Faint Source 
Survey (FSS) sources, Rowan-Robinson et al discov- 
ered IRAS F10214+4724, an IRAS galaxy with z = 
2.286 and a far infrared luminosity of 3x10 14 ^5 2 I/q . 
This object appeared to presage an entirely new class 
of infrared galaxies. The detection of a huge mass of 
CO by Brown and vandenBout (1991), 10 11 /i 50 2 M Q , 
confirmed by the detection of a wealth of molecular 
lines (Solomon et al 1992), and of submillimetre emis- 
sion at wavelengths 450-1250 /im (Rowan-Robinson 
et al 1991,1993, Downes et al 1992), implying a huge 
mass of dust, 10 h§ Mq confirmed that this was an 
exceptional object. Early models suggested this might 
be a giant elliptical galaxy in the process of forma- 
tion (Elbaz et al 92). Simultaneously with the grow- 
ing evidence for an exceptional starburst in F10214, 
the Seyfert 2 nature of the emission line spectrum 
(Rowan-Robinson et al 1991, Elston et al 1994a) was 
supported by the evidence for very strong optical po- 
larisation (Lawrence et al 93, Elston et al 94b). Subse- 
quently it has become clear that F10214 is a gravita- 
tionally lensed system (Graham and Liu 1995, Broad- 
hurst and Lehar 1995, Serjeant et al 1995, Eisenhardt 
et al 1996) with a magnification of about 10 at far 
infrared wavelengths, but not much greater than that 
(Green and Rowan- Robinson 1996). Even when the 
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magnification of 10 is allowed for, F10214 is still an 
exceptionally luminous far ir source. 

In 1992 Barvainis et al successfully detected sub- 
millimetre emission from the z=2.546 'clover-leaf 
gravitationally lensed QSO, H1413+117, which sug- 
gested that H1143 is of similar luminosity to F10214. 
Subsequently (Barvainis et al 1995) they realized that 
the galaxy was an IRAS FSC source. 

Here I want to place emphasis on the hyperlumi- 
nous galaxies detected as a result of unbiassed sur- 
veys at far infrared (and submillimctre) wavelengths. 
The program of follow-up of IRAS FSS sources which 
led to the discovery of F10214 (Rowan-Robinson et 
al 1991, Oliver et al 1996) has also resulted in the 
discovery of a further 6 galaxies or quasars with far 
ir luminosities > 10 13 ' 22 /i^ 2 L Q (McMahon et al 99). 
Four galaxies from the PSCz survey (Saunders et al 
1996) of IRAS galaxies brighter than 0.6 Jy at 60 /im 
fall into the hyperluminous category (a further two 
are blazars, 3C345 and 3C446, and these are not con- 
sidered further here), and a further example has been 
found by Stanford et al (1999) in a survey based on 
a comparison of the IRAS FSS with the VLA FIRST 
radio survey. Two galaxies detected in submillimetre 
surveys at 850 /im with SCUBA also fall into the hy- 
perluminous category (but one of these only because 
of the effect of gravitational lensing). 

Cutri et al (1994) reported a search for IRAS 
FSS galaxies with 'warm' 25/60 /jm colours, which 
yielded the z = 0.93 Seyfert 2 galaxy, F15307+3252 
(see also Hines et al 1995). Wilman et al (1999) have 
reported a further 2 hyperluminous galaxies from a 
more recent search for 'warm' galaxies by Cutri et al 
(1999, in preparation). Dey and van Breugel (1995) 
reported a comparison of the Texas radio survey with 
the IRAS FSS catalogue, which resulted in 5 galaxies 
with fir ir luminosities > 10 13 fe^" 2 Lq. However three 
of these are present only in the FSS Reject Catalogue 
and have not been able to be confirmed as far in- 
frared sources to date. The other two are discussed 
below. Four PC quasars from the list of Sanders et al 
(1989) (two of which are part of the study of Rowan- 
Robinson (1995)), fall into the hyperluminous cate- 
gory. Recently Irwin et al (1999) have found a z = 
3.91 quasar which is associated with with IRAS FSS 
source F08279+5255, the highest redshift IRAS ob- 
ject to date. 

Finally, inspired by the success in finding highly 
redshifted submillimetre continuum and molecular 
line emission in F10214, several groups have studied 
an ad hoc selection of very high redshift quasars and 
radio-galaxies, with several notable successes (An- 
dreani et al 1993, Dunlop et al 1995, Isaak et al 1995, 
McMahon et al 1995b, Ojik et al 1995, Ivison 1995, 
Omont et al 1997, Hughes et al 1997, McMahon et al 
1999). Many of these detections imply far ir luminosi- 
ties > lO 13 ' 22 /i5 O 2 Z/0 , assuming that the far ir spectra 
are typical starbursts. In all there are now 39 hyper- 
luminous infrared galaxies known, which are listed in 



Tables 1-3 according to whether they are (1) the re- 
sult of unbiased 60 fim (or submm) surveys, (2) found 
from comparison of known quasar and radio-galaxy 
lists with 60 /im catalogues, (3) found through sub- 
millimetre observations of known high redshift AGN. 
Table 4 lists some luminous infrared galaxies which 
do not quite meet my criteria, but have far infrared 
luminosities > W 130 Lq. But to set these in perspec- 
tive there are a further 20 PSCz galaxies which have 
13.00 < logwL ir /L Q < 13.22 (for H a = 50). 

From the surveys summarised in Table 1 we can 
estimate that the number of hyperluminous galaxies 
per sq deg brighter than 200 mjy at 60 /im is 0.0027- 
0.0043, which would imply that there are 100-200 hy- 
perluminous IRAS galaxies over the whole sky, 25 of 
which are listed in Tables 1 and 2. 



4 MODELS FOR HYPERLUMINOUS 
INFRARED GALAXIES 

For a small number of these galaxies we have rea- 
sonably detailed continuum spectra from radio to uv 
wavelengths. Figures 1-17 show the infrared continua 
of these hyperluminous galaxies, with fits using radia- 
tive transfer models (specifically the standard M82- 
like starburst model and an Arp220-like high optical 
depth starburst model from Efstathiou et al (1999) 
and the standard AGN dust torus model of Rowan- 
Robinson (1995). More than half of those shown have 
measurements at at least 9 independent wavelengths 

We now discuss the individual objects (and a few 
not plotted) in turn: 

F10214 + 4724 

The continuum emission from F10214 was the 
subject of a detailed discussion by Rowan-Robinson 
et al (1993). Green and Rowan-Robinson (1995) have 
discussed starburst and AGN dust tori models for 
F10214. Fig 1 shows M82-like and Arp 220-like star- 
burst models fitted to the sumbillimetre data for this 
galaxy. The former gives a good fit to the latest data. 
The 60 fj,m flux requires an AGN dust torus compo- 
nent. To accomodate the upper limits at 10 and 20 
/im, it is necessary to modify the Rowan-Robinson 
(1995) AGN dust torus model so that the maximum 
temperature of the dust is 1000 K rather than 1600 K. 
I have also shown the effect of allowing the dust torus 
to extend a further factor 3.5 in radius. This still does 
not account for the amplitude of the submm emission. 
The implied extent of the narrow-line region for this 
extended AGN dust torus model, which we use for 
several other objects, would be 326 (Lboi/W 13 Lq) 1 ^ 2 
pc consistent with 60-600 (Ltoi / W 13 Lq) 1 ^ 2 pc quoted 
by Netzer (1990). Evidence for a strong starburst con- 
tribution to the ir emission from F10214 is given by 
Kroker et al (1996) and is supported by the high gas 
mass detected via CO lines (see section 6) . Granato et 
al (1996) attempt to model the whole sed of F10214 
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with an AGN dust torus model, but still do not ap- 
pear to be able to account for the 60 /im emission. 
F0023 + 1024 

A starburst model fits the IRAS and ISO data 
(Verma et al 1999) well and there is a strong limit on 
any AGN dust torus contribution . 

SMMJ02399 0136 

A starburst model fits the submm data well and 
the ISO detection at 15 /im gives a very severe con- 
straint on any AGN dust torus component. The star- 
burst interpretation of the submm emission is sup- 
ported by the gas mass estimated from CO detections 
(Frayer et al 1999b). 

F1421 + 3845 

A starburst model is the most likely interpreta- 
tion of the IRAS and ISO 60-180 fim data, but there 
are discrepancies. There is a strong limit on any AGN 
dust torus component . 

TX0052 + 4710 

There is little evidence for a starburst contribu- 
tion to the sed of this galaxy. An extended AGN dust 
torus model fits the data reasonably well, apart from 
the ISO detection by Verma et al (1999) at 180 /im. 

F08279 + 5255 

An M82-like starburst is a good fit to the submm 
data and an AGN dust torus model is a good fit to 
the 12-100 /im data. The high gas mass detected via 
CO lines (Downes et al 1998) supports a starburst 
interpretation, though the ratio of L a b/M ga s is on the 
high side (see section 6). The submm data can also 
be modelled by an extension of the outer radius of 
the AGN dust torus and in this case the starburst 
luminosities given in Table 2 will be upper limits. 

P09104 + 4109 

The 100 /im upper limit places a limit on the 
starburst contribution and there is no detection of 
CO emission in this galaxy. The 12-60 /im data can 
be fitted by the extended AGN dust torus model. 

PG1148 + 549 

A starburst is a natural explanation of the 100 
/xm excess compared to the AGN dust torus fit to 
the 25 and 60 /im data, but detection in the submm 
and in CO would be important for confirming this 
interpretation. 

PG1206 + 459 

The IRAS 12-60 /tm data and the ISO 12-200 /im 
data (Haas et al 1998) are well-fitted by the extended 
AGN dust torus model and there is no evidence for a 
starburst. 

H1413 + 117 

The submm data is well fitted by an M82-like 
starburst and the gas mass implied by the CO detec- 
tions (Barvainis et al 1994) supports this interpreta- 
tion. The extended AGN dust torus model discussed 
above does not account for the submm emission. How- 
ever Granato et al (1996) model the whole sed of 
H1413 in terms of an AGN dust torus model. 

15307 + 325 

A starburst model gives a natural explanation 



for the 60-180 /im excess compared to the AGN dust 
torus model required for the 6.7 and 14.3 /tm emission 
(Verma et al 1999), but the non-detection of CO poses 
a problem for a starburst intepretation. 
PG1634 + 706 

The IRAS 12-100 /im data and the ISO 150-200 
/im data (Haas et al 1998) are well-fitted by the ex- 
tended AGN dust torus model and an upper limit 
can be placed on any starburst component. The non- 
detection of CO is consistent with this upper limit. 

4C0647 + 4134 

An M82-like starburst model gives a reasonable 
fit to the submm data (although the 1250 /im flux 
seems very weak). Observations in the mid-ir are 
needed to constrain the AGN dust torus. Since the 
AGN is not seen directly we have no constraints on 
its optical luminosity. We have used the non-detection 
of this source by IRAS (which we take to imply S(60) 
< 250 mjy) to set a limit on Ltor- This limit is not 
strong enough to rule out an AGN dust torus inter- 
pretation of the submm data. 

BR1202 0725 

The submm data are well-fitted by an M82-like 
starburst model. The QSO is seen directly, although 
with strong self-absorption in the lines (Storrie- Lom- 
bardi et al 1996), so we can probably not use the QSO 
bolometric luminosity to set a limit on Ltor- Wilkes 
et al (1998) have detected this galaxy at 7, 12 and 25 
/im with ISO, which would imply that the QSO is un- 
dergoing very strong extinction. An AGN dust torus 
model is capable of accounting for the whole spec- 
trum but the starburst interpretation of the submm 
emission is supported by the gas mass estimated from 
CO detections (Ohta et al 1996, Omont et al 1996). 

PG1241 + 176 

The ISO data of Haas et al (1999) can be fitted 
with an AGN dust torus model and there is no evi- 
dence for a starburst component. The 1.3 mm flux is 
probably an extrapolation of the radio continuum. 

PG1247 + 267 The ISO data of Haas et al 
(1999) can be fitted with an AGN dust torus model 
and there is no evidence for a starburst component. 

PG1254 + 047 The ISO data of Haas et al 
(1999) can be fitted with an AGN dust torus model 
and there is only weak evidence for a starburst com- 
ponent. 

BRI1335 0417 

The submm data can be fitted with a starburst 
model and since the QSO is seen directly we can use 
its estimated bolometric luminosity to set a limit on 
Ltor, which makes it unlikely that the submm emis- 
sion is from an dust torus. The starburst interpre- 
tation is supported by gas mass estimated from CO 
detections (Guilloteau et al 1997). 

PC2047 + 0123 

The limit on Ltor from the bolometric luminosity 
of the QSO makes it unlikely that an AGN dust torus 
is responsible for the 350 /im emission. However a 
starburst model can not fit both the 350 and 1250 /im 
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observed fluxes. Observations at other submm wave- 
lengths may help to clarify the situation. 

For the remaining objects we have only 60 fim or 
single submillimetre detections and for these we esti- 
mate their far infrared luminosity, and other proper- 
ties, using the standard starburst model of Efstathiou 
et al (1999). Tables 1-4 give the luminosities inferred 
in the starburst (L s b) (and fits of the A220 model 
in brackets) and AGN dust torus (L tor ) components, 
or limits on these, with an indication, from the row 
position of the estimate, of which wavelength the es- 
timate is made at. In Fig 18 we show the far infrared 
luminosity derived for an assumed starburst compo- 
nent, versus redshift, for hyperluminous galaxies, with 
lines indicating observational constraints at 60, 800 
and 1250 /im. Three of the sources with (uncorrected) 
total bolometric luminosities above 10 14 /i5 2 L o are 
strongly gravitationally lensed. IRAS F10214+4724 
was found to be lensed with a magnification which 
ranges from 100 at optical wavelengths to 10 at far in- 
frared wavelengths (Eisenhardt et al 1996, Green and 
Rowan- Robinson 1996). The 'clover- leaf ' lensed sys- 
tem H1413+117 has been found to have a magnifica- 
tion of 10 (Yun et al 1997). Downes et al (1998) report 
a magnification of 14 for F08279+5255. Also, Ivison 
et al estimate a magnification of 2.5 for SMMJ02399- 
0136 and Frayer et al (1999a) quote a magnification 
of 2.75 ±0.25 for SMMJ14011+0252. 

These magnifications have to be corrected for in 
estimating luminosities (and dust and gas masses) 
and these corrections are indicated in Fig 18. It ap- 
pears to be a reasonable assumption that if a star- 
burst luminosity in excess of 10 14 L Q is measured 
then the source is likely to be lensed, so F14218 and 
TX1011 merit further more careful study. 

On the other hand there is strong evidence for 
a population of galaxies with far ir luminosities in 
the range l-3xlO 13 /i^ o 2 I/0 . I have argued that in 
most cases the rest-frame radiation longward of 50 
/j,m comes from a starburst component. The luminosi- 
ties are such as to require star formation rates in the 
range 3-lOxlO 3 /i^ o 2 M0yr, which would in turn gener- 
ate most of the heavy elements in a 10 1 Mq galaxy 
in 10 7 — 10 s yrs. Most of these galaxies can therefore 
be considered to be undergoing their most significant 
episode of star formation, ie to be in the process of 
'formation'. 




Figure 1. Observed spectral energy distribution for 
F10214, modelled with M82-typc starburst (solid curve), 
Arp 220- type starburst (broken curve), AGN dust torus 
(dotted curve, modified AGN dust torus model - long- 
dashed curve). 




log 10 X(yum) 



5 THE ROLE OF AGN 

It appears to be significant that a large fraction of 
these objects are Seyferts, radio-galaxies or QSOs. For 
the galaxies in Tables 2 and 3, this is a selection ef- 
fect in that these objects are deliberately selected to 
be, or to be biased towards, AGN. For the popula- 
tion of objects found from direct optical follow-up of 
IRAS samples or 850 /itm surveys (Table 1), out of 
12 objects, 5 are QSOs or Seyfert 1, 1 is Seyfert 2, 



Figure 2. Observed spectral energy distribution for 
F0023, notation as for Fig 1. 



and 6 are narrow-line objects. Thus in at least 50 % 
of cases, this phase of exceptionally high far ir lu- 
minosity is accompanied by AGN activity at optical 
and uv wavelengths. This proportion might increase 
if high resolution spectroscopy were available for all 
the galaxies. For comparison the proportion of ultra- 
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SMMJ02399, z=2.8 
QSO 





Figure 3. Observed spectral energy distribution for 
SMMJ02399, notation as for Fig 1. 



Figure 5. Observed spectral energy distribution for 
TX0052, notation as for Fig 1. 
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Figure 4. Observed spectral energy distribution for 
F1421, notation as for Fig 1. 



Figure 6. Observed spectral energy distribution for 
F08279, notation as for Fig 1. 



luminous galaxies which contain AGN has also been 
estimated as 49 % (Veilleux et al 1999). However de- 
spite the high proportion of ultraluminous and hyper- 
luminous galaxies which contain AGN, this does not 
prove that an AGN is the source of the rest-frame far 
infrared radiation. The ISO-LWS mid-infrared spec- 
troscopic programme of Genzel et al (1998), Lutz et 
al (1998), has shown that the far infrared radiation 
of most ultraluminous galaxies is powered by a star- 



burst, despite the presence of an AGN in many cases. 
Wilman et al (1999) have shown that the X-ray emis- 
sion from several hyperluminous galaxies is too weak 
for them to be powered by a typical AGN. 

In the Sanders et al (1989) picture, the far in- 
frared and submillimetre emission would simply come 
from the outer regions of a warped disk surround- 
ing the AGN. Some weaknesses of this picture as 
an explanation of the far infrared emission from PG 
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Figure 7. Observed spectral energy distribution for 
P019104, notation as for Fig 1. 



Figure 9. Observed spectral energy distribution for 
PG1206, notation as for Fig 1. 





Figure 8. Observed spectral energy distribution for 
PG1148, notation as for Fig 1. 



Figure 10. Observed spectral energy distribution for 
H1413, notation as for Fig 1. 



quasars have been highlighted by Rowan-Robinson 
(1995). A picture in which both a strong starburst 
and the AGN activity are triggered by the same in- 
teraction or merger event is far more likely to be capa- 
ble of understanding all phenomena (cf Yamada 1994, 
Taniguchi et al 1999). 

Where hyperluminous galaxies are detected at 
rest- frame wavelengths in the range 3-30 fim (and this 
can correspond to observed wavelengths up to 150 



/im), the infrared spectrum is often found to corre- 
spond well to emission from a dust torus surrounding 
an AGN (eg Figs 1, 6-10, 12). This emission often con- 
tributes a substantial fraction of the total infrared (1- 
1000 /im) bolometric luminosity. For the 12 ir-selected 
objects of Table 1, the luminosity in the dust torus 
component exceeds that in the starburst for 5 of the 
galaxies (42%). The advocacy of this paper for lu- 
minous starbursts relate only to the rest-frame emis- 
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Figure 11. Observed spectral energy distribution for 
F15307, notation as for Fig 1. 



Figure 13. Observed spectral energy distribution for 
4C0647, notation as for Fig 1. 
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Figure 12. Observed spectral energy distribution for 
PG1634, notation as for Fig 1. 



Figure 14. Observed spectral energy distribution for 
BR1202, notation as for Fig 1. 



sion at wavelengths > 50pim. Figure 19 shows the 
correlation between the luminosity in the starburst 
component, L a b, and the AGN dust torus component, 
Ltor, for hyperluminous infrared galaxies, PG quasars 
(Rowan- Robinson 1995), and IRAS galaxies detected 
in all 4 bands (Rowan-Robinson and Crawford 1989) 
(this extends Fig 8 of Rowan- Robinson 1995). The 
range of the ratio between these quantities, with 0.1 
< L 3 b/Ltor < 10, is similar over a very wide range of 



infrared luminosity (5 orders of magnitude), showing 
that the proposed separation into these two compo- 
nents for hyperluminous ir galaxies is not at all im- 
plausible. 



6 DUST AND GAS MASSES 

The radiative transfer models can be used to derive 
dust masses and hence, via an assumed gas-to-dust ra- 
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Figure 18. Bolometric luminosity in starburst component for galaxies with luminosities > 10 13 Lq (Tables 1-4). The 
galaxies labelled L are known to be lensed. Loci corresponding to the limits set by S(60) = 200 mjy, S(800) = 20 mjy, 
and S(1250) = 2 mjy are shown. 



tio, gas masses. For the M82-like starburst model used 
here the appropriate conversion is Mdust = 1CP 4 : ' 6 L s b, 
in solar units (Green and Rowan- Robinson 1996). 
These estimates have been converted into estimates 
of gas mass assuming M gaB jMduat = 300 (tables 1- 
4, col 9, bracketed values). However these estimates 
will not assist in deciding the plausibility of the star- 
burst models, because the radiative transfer mod- 
els are automatically self-consistent models of mas- 
sive star-forming molecular clouds. Estimates derived 
from v^BviTd) fits to the spectral energy distribu- 
tions are even less physically illuminating. 

Far more valuable are the cases where direct esti- 
mates of gas mass can be derived from molecular line 
(generally CO) observations. Where available, these 
estimates have been given in Tables 1-4, col 9, taken 
from Frayer et al 1999a (and references therein), Bar- 
vainis et al 1998, Evans et al (1999), Yun and Scoville 
(1999). Figure 20 shows a comparison of the estimates 
of L sb derived here with estimates of M gaa derived 
from CO observations. Also shown are results for ul- 
traluminous ir galaxies (Solomon et al 1997) and for 



more typical IRAS galaxies (Sanders et al 1991) (after 
rationalisation of some objects in common). 

The appropriate conversion factor from CO lu- 
minosity to gas mass is a matter of some controversy. 
For luminous infrared galaxies, Sanders et al (1991) 
used a characteristic value for molecular clouds in 
our Galaxy, 4.78 M e (Kkms~ 1 pc~ 2 y 1 . Solomon et 
al (1997) found that such a value led to gas mass es- 
timates for ultraluminous infrared galaxies a factor 
of 3 or more in excess of the dynamical masses and 
concluded that a value of 1.4 M a (Kkms~ 1 pc~ 2 )~ 1 
was more appropriate for these galaxies. Downes 
and Solomon (1998) studied several ultraluminous in- 
frared galaxies in detail in CO 2-1 and 1-0 with the 
IRAM interferometer, deriving an even lower value 
of 0.8 MQ(Kkms~ 1 pc~ 2 )~ 1 on the basis of radia- 
tive transfer models for the CO lines. However their 
gas masses are on average only l/6th of the (re- 
vised) inferred dynamical masses. In their detailed 
model for Arp 220, Scoville et al (1997) found a con- 
version factor 0.45 times the Galactic value, ie 2.15 
Mo(Kkms~ 1 pc~ 2 )~ 1 . Combining this with an esti- 
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Figure 19. Bolomctric luminosity in AGN dust torus component versus bolometric luminosity in star burst component: 
filled circles, hyperluminous ir galaxies (this paper); open triangles, PG quasars (Rowan-Robinson 1995); crosses, IRAS 
galaxies detected in 4 bands (Rowan- Robinson and Crawford 1989, galaxies with only upper limits on Ltor omitted). 



mated ratio for T(3-2)/T(l-0) of 0.6, Frayer et al 
(1999a) justify a value of 4 M e (Kkms~ 1 pc~ 2 )~ 1 for 
gas mass estimates of hyperluminous galaxies derived 
from CO 3-2 observations. 

In Tables 1-4 and Fig 19, I have followed 
Frayer et al (1999a) in using a conversion factor 
of 4 M El (Kkms~ 1 pc~ 2 )~ 1 for hyperluminous galax- 
ies. For other galaxies in Fig 19 with luminosities 
> 10 n ' 5 L Q , I have used a conversion factor of 2 
M Q {Kkms' 1 pc' 2 y 1 (in line with Scoville et al 1997, 
but a factor of 2 or so higher than advocated by 
Downes and Solomon 1998); for galaxies with lumi- 
nosities < 1O 11,5 L , I have used the standard Galac- 
tic value, 4.78 M Q (Kkms~ 1 pc~ 2 )~ 1 . 

The range of ratios of L s b/M gas for hyperlumi- 
nous galaxies is consistent with that derived for ul- 
traluminous starbursts. For cases where we have esti- 
mates of gas mass both from CO lines and from dust 
emission, the agreement is remarkably good (within 
a factor of 2). There is a tendency for the time-scale 
for gas-consumption, assuming a star formation rate 
given by eqn (1), to be shorter for the more luminous 



objects, in the range 10 7 — 10 8 yrs (alternatively this 
could indicate a higher value for the low-mass cutoff 
in the IMF). The cases where a strong limit can be set 
on Mgas are also, generally, those where the seds do 
not support the presence of a starburst component. 
After correction for the effects of gravitational lens- 
ing, gas masses ranging up to 1-3 x 10 11 Mq are seen 
in most hyperluminous galaxies, comparable with the 
total stellar mass of an L* galaxy (10 112 (M/4L)/i5 2 ). 
In fact 24/39 hyperluminous galaxies in Tables 1-3 
have gas masses estimated either from CO or from 
dust emission > 1O 11 M (after correction for effects 
of lensing, where known). Hughes et al (1997) ar- 
gue that a star-forming galaxy can not be consid- 
ered primeval unless it contains a total gas mass of 
1O 12 M0, but this seems to neglect the fact that 90 % 
of the mass of galaxies resides in the dark (probably 
non-baryonic) halo. 
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Figure 20. Bolometric luminosity in starburst component, versus mass in gas, deduced from CO observations. Filled 
circles, hypcrluminous ir galaxies (Frayer et al 1999, Yun and Scovile 1999); open triangles, ultraluminous ir galaxies 
(Solomon ct al 1997); crosses, IRAS galaxies (Sanders ct al 1991). 



7 CONCLUSIONS 

(1) About 50 % of hyperluminous infrared galaxies se- 
lected in unbiassed infrared surveys have AGN optical 
spectra. This is not in fact higher than the propor- 
tion seen in ultraluminous ir galaxies by Vielleux et 
al (199a). For about half of the galaxies in this sample, 
the AGN dust torus is the dominant contribution to 
the total ir (1-1000 fim) bolometric luminosity, while 
in half of cases a starburst seems to be the dominant 
contributor. 

(2) There is a need for both an AGN dust torus 
and starburst components to understand most seds 
of hyperluminous ir galaxies (29/39). Measured gas 
masses support, in most cases, the starburst interpre- 
tation of rest-frame far-infrared and submm (A em > 
50^im)emission. 

(3) There is a broad correlation between the lu- 
minosities of starburst and AGN dust torus compo- 
nents (Fig 19). This may imply that there is a physical 
link between the triggering of star formation and the 
feeding of a massive black hole. Taniguchi et al (1999) 
have argued that during a merger giving rise to a lu- 



minous starburst, a pre-existing black hole of 10 M© 
may grow into a large one > 1O 8 M0 and hence form a 
quasar. Alternatively they suggest that a large black 
hole might be formed out of star clusters with com- 
pact remnants. 

(3) There is no evidence in most objects that an 
AGN powers a significant fraction of radiation at rest- 
frame wavelengths > 50 urn. For P09104 and PG1634, 
the non-detection of CO emission is consistent with 
the absence of evidence in the sed for a starburst 
component. In F08279, the mass of CO detected sug- 
gests a limit on the starburst luminosity which implies 
that the observed submm radiation may simply be the 
long-wavelength tail of its AGN dust torus emission. 
F15307 poses a problem: the sed can be understood 
as radiation from both an AGN dust torus and an 
Arp-220 like starburst, but the upper limit on the 
molecular mass from the non-detection of CO would 
then imply a very extreme ir-luminosity to gas-mass 
ratio. 

(4) After correction for the effects of lensing, star- 
formation rates in excess of 2000 Mq /yr are inferred 
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Figure 15. Observed spectral energy distribution for 
PG1241+176, notation as for Fig 1. 



Figure 17. Observed spectral energy distribution for 
PG1254+047, notation as for Fig 1. 
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Table 1. Hypcrluminous Infrared Galaxies, found in 60 and 850 fim surveys 



P00182-7112 



F0023+1024 



P07380-2342 



F10026+4949 



F10214+4724 



F12509+3122 



F1327+3401 



P 14026+4341 



F14218+3845 



0.327 



0.58 



SMMJ02399-0136 2.803 



0.292 



1.12 



2.286 



0.78 



0.36 



0.324 



1.21 



A 


fiux(mjy) 


ref 


sp.type lg L sb lg L tor 


100 


1120. 


32 


nl 


60 


1300. 




13.25 (13.32) 


25 


<250. 




<13.06 


12 


<250. 






180 


1047 ±52 


1,23 


nl 


100 


<938 






90 


923 ±78 


23 




60 


428 




13.37 (13.34) 


25 


<193 




< 12.23 


14.3 


<10.2 


23 




12 


<173 






6.7 


<2.45 


23 




1350 


5.7 ±1.0 


25 


S2(L) 


850 


26 ± 3 




13.52 (13.11) 


750 


28 ±5 






450 


69 ±15 






350 


<323. 






100 


<715. 






60 


<428. 






25 


<86. 






15 


1.2 ±0.4 




13.10 


12 


<91. 






100 


3550. 


32 


nl 


60 


1170. 




13.09 (13.15) 


25 


800. 




13.46 


12 


480. 






100 


<619 


1 


SI 


60 


266 




<13.84 (<13.90) 


25 


(177) 




14.04 


12 


<86 






1200 


9.6 ± 1.4 


11 


S2(L) 


1100 


24 ± 5 


7 




800 


50 ± 5 




13.83 (13.45) 


450 


273 ± 45 






350 


383 ±51 


29 




100 


<510 






60 


190 ± 40 




14.41 


20 


<45 






10 


<12 






100 


<675 


1 


QSO 


60 


218 




13.43 (13.63) 


25 


(103) 




13.51 


12 


<106 






100 


(1200) 


1 


QSO 


60 


1180 




13.30 (13.39) 


25 


<126 




< 12.85 


12 


<94 






100 


(1280.) 


32 


QSO 


60 


610. 




12.91 (12.99) 


25 


260. 




13.07 


12 


<450. 






180 


1575 ±115 


1,23 


QSO 


100 


(2100) 






90 


1985 ±129 


23 




60 


565 


1 


14.39 (14.71) 


60 


1335 ±255 


23 




25 


<75 






14.3 


<5.4 


23 




12 


<97 






6.7 


1.07 ±0.37 


23 


<12.65 



(11.13) 



(11.25) 



11.65 (11.40) 



(10.97) 



11.63 (11.71) 



(11.31) 



(11.18) 



(10.79) 



(12.27) 
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Table 2. Hyperluminous Infrared Galaxies, found in 60 and 850 fim surveys (cont.) 



name 


z 


A 


flux(mjy) 


ref 


sp.type 




lg Ltor 


lg Mgas 


P18216+6418 


0.30 


100 


2160. 


32 


nl 












60 


1130. 






13.10 (13.17) 




(10.98) 






25 


400. 
















12 


(190.) 








13.15 




FFJ1614+3234 


0.710 


100 


<540. 


33 


nl 












60 


174. 






13.22 (13.40) 




(11.10) 






25 


<55. 








<13.15 








12 


<65. 












F2356-0341 


0.59 


180 


<707 


1,23 


nl 












100 


<792 
















90 


<251 


23 














60 


347 






13.30 (13.45) 




(11.18) 






25 


<142 
















14.3 


<9.3 


23 














12 


<87 
















6.7 


<1.74 


23 






<12.08 
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Table 3. Hypcrluminous Infrared Galaxies, found by comparison of 60 /im surveys with quasar or radio-galaxy lists, or 
using an infrared colour selection biased to AGN 



TX0052+4710 1.93 



F08279+5255 3.91 



P09104+4109 0.44 



TX1011+1430 1.55 



PG1148+549 



0.969 



PG1206+459 1.158 



PG1248+401 1.030 



H1413+117 



2.546 



F14481+4454 



0.66 



F14537+1950 0.64 



A 


flux(mjy) 


ref 


850 


<9.8 


2,22 


450. 


<128 


22 


180 


350 ±61 


23 


90 


277 ±84 


23 


60 


160 


2 


25 


<66 


2 


14.3 


<9.6 


23 


12 


<78 


2 


6.7 


<2.49 


23 


1350 


24 ±2 


26 


850 


75 ± 4 




450 


211 ±47 




100 


(951) 




60 


511 




25 


(226) 




12 


<101 




100 


<438 


6 


60 


525 




25 


333 




12 


(130) 




100 


<757 


2 


60 


225 




25 


<269 




12 


<116 




100 


410 


9 


60 


196 


9 


25 


120 


9 


12 


<75 


9 


170 


188 


36 


100 


386 


36 


60 


463 


36 


25 


<113 


9 


20 


208 


36 


12.8 


110 


36 


12 


207 ±36 


9 


4.8 


<10 


36 


100 


<378 


9 


60 


224 ±51 


9 


25 


< 200 


9 


12 


<U7 


9 


1250 


18 ± 2 


21 


800 


66 ± 7 


18 


761 


44 ± 8 


12 


438 


224 ± 38 


12 


345 


189 ± 56 


12 


350 


293 ±14 


29 


200 


280 ±50 


36 


100 


370 ± 78 


21 


60 


230 ± 38 


21 


25 


75 ±3 


36 


12 


57 ±8 


36 


7 


17 ±2 


36 


4.8 


10 ±3 


36 


100 


<500 


27 


60 


190 




25 


(85) 




12 


<76 




100 


<738 


27 


60 


283 




25 


< 159 




12 


< 119 





sp.typc 



<13.15 (<12.73) 



lg L to 



lg Mgas 

(<11.03) 



14.17 



<13.93 (<13.61) 



11.55 (<11.81) 



15.30 

S2 

<12.78 (<12.70) <10.48 (<10.66) 

13.39 

QSO 

14.32 (14.72) < 14.13 (12.20) 

QSO 13.50 (13.53) (11.38) 

13.74 

QSO 



<13.13 (<13.00) 



(<11.01) 



14.21 



QSO <13.53 (<13.57) (<11.41) 

13.83 



BALQ(L) 



13.94 (13.54) 



12.07 (11.82) 



14.61 



S2 



sb 



13.17 (13.34) (11.05) 
13.26 



13.31 (13.46) (11.19) 
< 13.51 
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Table 4. Hyperluminous Infrared Galaxies, found by comparison of 60 fim surveys with quasar or radio-galaxy lists, 
using an infrared colour selection biased to AGN 

sp.type lg L sb lg Ltor lg M gas 

F15307+325 0.93 3000. <1.3 31 S2 

13.55 (13.12) <10.39 (11.43) 

13.08 



A 


flux(mjy) 


ref 


3000. 


<1.3 


31 


1250 


<5.1 


31 


180 


397 ±67 


23 


100 


510 ± 62 


16 


90 


478 ±65 


23 


60 


280 ± 42 




25 


80 ± 24 




14.3 


15.3 ±2.6 


23 


12 


<45 




6.7 


2.95 ±0.51 


23 


800 


<47 


17 


450 


<420 


17 


200 


216 


35 


170 


238 


35 


150 


307 


35 


100 


343 


9 


60 


318 


9 


25 


147 


9 


12 


61 


9 



PG1634+706 1.334 800 <47 17 QSO 



<13.20 (<13.05) <10.70 (<11.08) 

14.13 
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Table 5. Hypcrluminous Infrared Galaxies, found through submm observations of known high rcdshift AGN 



name 


z 


A 


flux(mjy) 


ref 


sp.type 




lg L tor 


lg M gas 


HM0000-263 


4.10 


350 


134 ±29 


29 


QSO 


13.87 (13.79) 


<15.07 


(11.75) 


Q0100+1300 


2.68 


350 


131 ±28 


29 


QSO 


13.60 (13.42) 


< 14.70 


(11.48) 


PC0307+0222 


4.379 


1250 


6.6 ± 1.7 


4 


QSO 


13.31 (12.92) 


<13.23 


(11.19) 






R 


20.39 


20 










PC0345+0135 


3.638 


1250 


6.1 ± 2.0 


4 


QSO 


13.34 (12.93) 


<13.26 


(11.22) 






800 


<25 


18 














R 


19.94 


20 










MG0414+0534 


2.639 


350 


<105. 


29,34 


QSO(L) 












60 


140. ±40. 






<14.76 (<15.3) 


14.43 


11.53 (<12.64) 






25 


70. ±24. 












4C0647+4134 


3.8 


1250 


2.5 ± 0.4 


5 


RG 








(4C41.17) 




800 


17.4 ± 3.1 


5 




13.30 (12.97) 


<15.00 


(11.18) 






450 


<56 


5 














350 


37 ±9 


29 










BR1202-0725 


4.69 


1250 


10.5 ± 1.5 


18 


QSO 


<13.67 (<13.39) 


15.50 


11.41 (<11.55) 






1100 


21 ± 5 


8 














850 


42 ± 2 


8 














450 


92 ± 38 


8 














350 


106 ± 7 


29 














25 


165 ±25 


38 














12 


105 ±7 


38 














7 


15 ±2 


38 














4 


<6 


38 














R 


18.7 


20 










PG1241+176 


1.273 


1300 


3.3 


37 


QSO 












170 


<96 






<12.92 (<12.80) 




(<10.80) 






100 


217 
















60 


378 








14.28 








12.8 


<87 
















7.3 


<27 
















4.8 


7 












PG1247+267 


2.038 


1300 


<2.1 


37 


QSO 


<13.13 (<12.74) 




(<11.01) 






170 


<150 
















100 


174 
















60 


236 








14.26 








12.8 


30 
















7.3 


10 
















4.8 


<9 












PG1254+047 


1.024 


1300 


<1.8 


37 


QSO 


<13.14 (<12.75) 




(<11.02) 






170 


<345 
















100 


242 
















60 


307 








14.28 








12.8 


<90 
















7.3 


33 
















4.8 


27 












LBQS1230+1627 


2.735 


1250 


7.5 ± 1.4 


19 


QSO 


13.54 (13.05) 


< 14.31 


(11.42) 






350 


104 ±21 


29 














V 


17.4 


30 










BRI1335-0417 


4.396 


1250 


10.3 ± 1.04 


19 


QSO 


13.50 (13.12) 


<13.63 


11.44 (11.38) 






850 


14 ± 1 


24 














350 


52 ± 8 


29 














R 


19.4 


20 










PC2047+0123 


3.799 


1250 


2.08 ± 0.47 


13 


QSO 


13.59 (13.49) 


<13.39 


(11.47) 






350 


80 ±20 


29 














R 


19.71 


30 










PC2132+0126 


3.194 


1250 


11.5 ± 1.7 


4 


QSO 


13.66 (13.19) 


<13.22 


(11.54) 






800 


<12 


18 














R 


19.78 


30 
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Table 6. Luminous Infrared Galaxies not meeting the requirements for Tables 1-3, but satisfying L s f,+Lt or > 
10 13.0 



name 


z 


A 


flux(mjy) 


ref 


sp.type 


Ig L sb 


lg Ltor 


lg M gas 


BR0902+34 


3.391 


1250 


3.1 ±0.6 


18 


RG 


13.07 


<14.92 


(10.95) 






850 


<14 
















450 


<99 












BRI0952-0115 


4.43 


1250 


2.78 ± 0.63 


19 


QSO(L) 












850 


14 ± 2 


24 




13.19 (12.91) 


<13.92 


(11.07) 






350 


<66 


29 














R 


18.7 


20 










BR1117-1329 


4.00 


1250 


4.09 ± 0.81 


19 


QSO 


13.16 (12.86) 


<14.11 


(11.04) 






850 


13 ± 1 


24 














350 


<39 


29 














R 


18.0 


20 










F1220+0939 


0.68 


60 


180 


1 


QSO 


13.18 (13.36) 


<13.59 


(11.06) 


F12514+1027 


0.30 


100 


755. 


27 


S2 












60 


712. 






12.90 (12.97) 




(10.78) 






25 


190. 








12.86 








12 


<63.2 












SMMJ14011+0252 


2.55 


1350 


6.06 ±1.46 


28 


nl(L) 












850 


14.6 ±1.8 






13.36 (12.97) 


< 14.64 


11.49 (11. 






450 


41.9 ±6.9 
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